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Abstract 



' We investigate charge accretion in vicinity of a slowly rotating compact 

' star with a non-stationary electromagnetic field. Exact solutions to the 

. general relativistic Maxwell equations are obtained for a star formed of a 

' highly degenerate plasma with a gravitational field given by the linearized 

fH ^ Kerr metric. These solutions are used to formulate and then to study 

numerically the equations of motion for a charged particle in star's vicinity 
using the gravitoelectromagnetic force law. The analysis shows that close 
to the star charge accretion does not always remain ordered. It is found 
that the magnetic field plays the dominant role in the onset of chaos near 
52 ' the star's surface. 



- - ■ 1 Introduction 



Keywords: Plasma accretion; compact stars; chaos; chaotic dynamics; mag 
netic field; relativistic MHD, gravitoelectromagnetic effects 



Accretion is a fundamental process by which stars accumulate matter from their 
surroundings, it therefore plays an important role in stellar formation and evo- 
lution. For normal stars accretion is principally due to the gravitational field of 
the star. However for compact objects (classed as neutron stars, pulsars, and 
white dwarfs) accretion is not only due to gravitational but also electromagnetic 
effects. In vicinity of these stars matter is mostly in the form of a highly con- 
ducting electron plasma, whereas the star itself is formed of a degenerate highly 
conducting quasi-neutral plasma. The compact state of the plasma inside the 
star is brought about by the intense gravitational field of the star. Exterior to 
the star gravitational and electromagnetic field produce combined effects on the 
dynamics of plasma surrounding the star. Here gravitational field has a two-fold 



1 



effect on charge accretion. Firstly it effects the charge dynamics through the 
interaction of particle mass with the gravitational field of the star. Secondly 
the background spacetime also influences the electromagnetic field of the star 
which in turn effects the particle orbit via charge interaction. In general, under 
such conditions, the electromagnetic field of a compact star is not stationary 
but gains in intensity as the star evolves. Observationally the magnetic field 
intensity for a newly formed compact star to the final stages of a compact star 
can vary from lO'^G to lO^^G and hence varies with time'^"^]. 

In this paper we investigate the combined effects of gravitational and a non- 
stationary electromagnetic field on a charge particle dynamics in vicinity of a 
compact star. The star is considered to be formed of a degenerate highly con- 
ducting plasma. The magnetohydrodynamics (MHD) approximation is applied 
for the plasma forming the star. Using a generalized definition!'*''''! of the elec- 
tromagnetic field tensor for an MHD plasma in a curved spacetime the Maxwell 
field equations are explicitly formulated for a non-stationary axially symmet- 
ric electromagnetic field. Here the linearized Kerr metric is used to describe 
the gravitational field outside the star with ZAMO (zero angular momentum 
observers) as the stationary observer 1^1. The generalized Maxwell equation are 
then solved for the electric and magnetic field components using the separation 
of variable technique. These solutions are then used to obtain the equations of 
motion (geodesic equations) in linearized Kerr spacetime for a charge particle in 
the non-stationary electromagnetic field of the star. The numerical solutions to 
these equations show that at large distances from the star particle trajectories 
are well defined and ordered, however as the particle falls towards the star its 
orbit cannot be precisely predicted hence exhibits chaos. The solution show that 
star's magnetic field plays the dominant role in bringing the charge accretion to 
a chaotic state. 

The layout of the paper is as follows. In the next section generalized Maxwell 
equations are formulated using the generalized electromagnetic field tensor. 
These equations are solved in section 3 for a non-stationary axially symmetric 
electromagnetic field using a separation of variable ansatz for the field compo- 
nents. Then in section 4 equations of motion are formulated and numerically 
solved for a charged particle lying in the field of the star. In the conclusion 
(section 5) numerical results of section 4 are interpreted and the main results 
of the study are summarized. Throughout we have used the gravitational units 
in which G = 1 = c unless mentioned otherwise. 

2 Formulation of the Electromagnetic Field Equa- 
tions 

General relativistically the gravitational field of a rotating star can be described 

by the Kerr metric. In the region exterior to a slowly rotating compact star of 
mass M a linearization of the Kerr metric is given by 
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ds'' = -e'^'^^'^Ut^ -2uj{ry smHdtd<p + r'' sinHdif'' + e-^'^^''Ur^ +r^de^, (1) 

where we have used the usual Boyer-Lindquist coordinates. Also in expression 

(1) e^*'^''-' and w(r) are given by the relations e^*^''^ = (1 — 2M/r),uj{r) = 
d(fi/dt = —gtipl Qifip, where Lo{r) is the angular velocity of a free falling frame 
brought into rotation by the frame dragging of the spacetimel''! . The four ve- 
locity components tir and tig of a stationary observer ZAMO circling the star 
at a fixed radial distance vanish, therefore using = — 1, the components 

of the four velocity vector for a ZAMO are: 

«° = e-*W(l,0,0,u;(r)), 7X„ = e*M (-1, 0, 0, 0). (2) 

For a star formed of highly degenerate plasma the electromagnetic field (i?", B") 
of the star can be determined by the Maxwell equations. In a curved spacetime, 
such as the Kerr spacetime the general relativistic form of the Maxwell equations 
is given by 

■Pa/3,7 + ■P/37,a + -p7a,/3 = 0, (3) 

{^F-f')^^ = An^J'^, (4) 

where g represents the determinant of the metric tensor gap. Here Fap is the 
generalized electromagnetic field tensor for an ideal magnetohydrodynamic fluid 
given by a unique tensor expression: 

Fap = UaEfj - UfjEa + TJafS^su'' , (5) 

and J" is current four vector. In general the current four vector is the sum of 
two terms corresponding to a convection and to a conduction current: J" = 
eu" + au/sFl^", where e is the proper charge density and a is the conductivity 
of the plasma. The volume element 4-form r]a0-ys is equal to ^/^Cafj-ys where 
^ap-yS is the Levi -Civita symbol. Since in a compact star the quasi-neutral 
degenerate plasma is essentially a perfect conductor, we assume that for the 
region external to the star J" can be approximated to zero (see references 
[8] and [9] for a discussion of this often made assumption for a quasi-neutral 
degenerate plasma). 

To determine the non-stationary electromagnetic field of the star we as- 
sume that the electric and magnetic field components are dependent not only 
on the polar coordinates r and 9 but also the time coordinate t in the frame 
of ZAMO. We then obtain from the Maxwell equations (3) and (4) the fol- 
lowing determining equations for the electromagnetic field in the slow rotation 
approximationl^^l : 

drir"^ smeu^B"-) + dgir'^u* sineB^) = 0, (6) 

sm9u*dtB'^ - driutE^) = 0, (7) 
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sin eu^dtB'' + UfdeE^ = 0, (8) 
sin eu^dtB"^ + driutEg - sin Ou'^B'-) - deiufEr - sin fi'u'^S^) = 0, (9) 
dr{r^ smOu^E^) + deir'^u* sinOE^) = 0, (10) 
sinOu^dtE^ - dr{utB^) = 0, (11) 
sin eu^dtE'' + utdeB^ = 0, (12) 
sin6'u*at£;^ + dr{utBe - smOu'^E'') - de{utBr - sinOu'^E'^) = 0. (13) 

3 Solution for the Electromagnetic Field Exte- 
rior to the Star 

The non-stationary electromagnetic field of the star is now determined by the 
system of equations (6) to (13). To obtain the time dependent electric and mag- 
netic field components E°' and B" we assume the following separation ansatz 
for the field components: 



E^=A{r)ait)^ E'=B{r)b{t)^ E^=C{r)c{t)^ (14) 
sm 6 smo sm 6 



S'-=P(r)£.(0^, B'>=S{r)c{t)^ = J^{r)m^ (15) 

sin sin U sin 

Substituting into (6) we obtain 

^{r^u'V)^0 (16) 
dr 

This implies that ^{r) = const. /r^u*^. Also from equation (10) we obtain A{r) = 
const. /r^u^ hence the radial components of the electric and magnetic field reduce 
as: 

E''{r,e,t)=5^{t)/r'^u'sm^9, B'' {r,0,t) = aa{t)/r'^u* sin^ (17) 
From (13) it therefore follows that 

r'^u^Cj^ + ^{r^UtSte - r^uju^Aaa) - Ut smO^^ = (18) 
at dr do 

Comparing coefficients of sin 9 on both sides the last equation implies that B^ 
must be independent of 9. This is compatible with equation (17) iff (5 = 0. 
Therefore it follows that B^ = 0. Similarly from (17) it follows that E'' = 0. 
We therefore have after the separation of variables, from equation (13), (11), 
(9), and (7) respectively, the following determining equations 

1 dc — £ c? , 9 ^, , „. 
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fdi = W^drr^^^'-^^='^ (20) 

Idf -13 d 



b dt Xr^u^J^ dr 
1 dt 7 d 



(uyB) = -/C3 (21) 



where ki, k^, ks, and ^4 are the separation constants. 

3.1 TIME DEPENDENT SOLUTIONS 

For the separated equations (19) to (22) we have for the time dependent parts 
b{t), c{t), e{t), and j{t) the following set of equations: 

Differentiating the first of equations (23), and substituting in the second equa- 
tion of the pair (24) we obtain 

c{t) = Ai e^p{-^^/hklt) (25) 
Substituting back in equation (23) we obtain for e: 



-i^/kik 
ki 

where Ai is a constant. In a similar manner we obtain for f{t) and b{t): 



z(t) = , '^ Ai eM-iy/hk^t) (26) 
ki 



b{t) = Ai e^p{-i^/k^t) (27) 



^ »V^ ^^ exp(-iv/fc^i) (28) 
k2 

3.2 RADIAL SOLUTIONS 

The radial components C(r), £{r), T{r) and B{r) arc obtainable explicitly from 
equations (19) to (22) as follows. First from the radial part of the separated 
equation (19) we have for the component C(r): 



= \ 2 T^^ \^ ^) (29) 

Ki7r^ dr V r 

Substituting into equation (22) gives an expression for £{r) 
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^[il-—)^iJl-—r^S)]-^^^ (30) 
dr r dr \ r 2m 



This differential equation can be solved explicitly for the component £{r) to 
give 

S{r) = T-^l'^ir - 2M)-5+2Mx/fc7fcIc'i exp(yfclfclr)- 

r-3/2(r-2M)-i-2^^/*^— ^=exp(-\/fc^r) (31) 

The first part of the above solution is asymptotically unbounded, hence we take 
only the second part of the solution for E(r). Putting then for E(r) in equation 
(29) we obtain after some simplifications the component C(r): 



C(r) = ^ 3^(r-2M)-2-2A^yfcTfcI{^(2+^^r)-3M}exp(-V^r) 

(32) 

In a similar manner we obtain for J^{r) and B{r): 

r-3/2(r-2M)-^-2^^^-^exp(-^fc^r) (33) 



B{r) = ^- ^(r-2M)-5-2^v^{r(2+Vfc^r)-3M}exp(-v/^r) 

2/3r9/2v^^3 

(34) 

These solutions are single valued and bounded for r — > oo. 

Hence the radial solutions (31) to (34) together with the time dependent 
solutions (25) to (28) determine the electromagnetic field of a compact star in 
the region exterior to the star. 

4 Charge Accretion in the Slow Rotation Ap- 
proximation 

A linearization of the geodesic equation for a slowly rotating gravitational source 
is given by the gravitoelectromagnetic force of magnitude G + v x H per unit 
mass where G and H arc the gravitoelectric (GE) and gravitomagnetic (GM) 
potentialsl^^'^^i. In the slow rotation limit the star can be considered as a 
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homogeneous sphere of mass M and radius R. In this case the GE and GM 
potentials are given by 



G = -^r, H = -fMi^^(0.r^-i^l), (35) 



where 12 is the angular velocity vector of the star. In the presence of an elec- 
tromagnetic field the force law for a particle with mass m and charge q is given 
by 



— = (G + vxH)+^(E + vxB). (36) 



— = G + vxH 

at'^ m, 

The equations of motion can be studied mimerically iii Cartesian coordinates 
by transforming the polar coordinates into the Cartesian coordinates as x = 
r sin 9 cos (p,y = r sin 9 sin (p,z = r cos 9. The transformation gives for the elec- 
tric field components = — rsini^i?*', Ey = r cos ipE'^, = —riim9E^, 
and similarly for the magnetic field i?^ = —r sin ipB'^ , = rcos</?_B'^, = 
—rsm9B^. In a two dimensional plane 9 = it/2 (the equatorial plane of the 
star), the gravitomagnetic effects are maximum, here we get for the equations 
of motion for the (a;, y) plane: 

If = _ ^. + („ _ = + + ,H + ,B.)| (37) 

§ = -» + - (H - = -» + .E^ - iH + (38) 

where from (14) and (15) we have for the real part of -E*^ and B^ 



/l _ 2M 

Re E'P = Eo ^ (r - 2M) - 5 -2^*= {r (2 + fcr ) - 3M} exp(- fcr) cos kt (39) 

ReB^ = Bor-"^^^(r-2M)-3-2Mfcexp(-fcr)sinfci (40) 

where r = y/ x'^ + y'^ , k = \/^i^4j and E'o = — £AiC2/2fcifc4 sin^ 6* and _Bo = 
eAiC^ /2fci sin are constants for fixed value of 9, such as in the equatorial plane 
of the star. Also g=\G\= M/{x^+y'^), and H =\'R\= n/{x^+y'^f/'^, whereas 
H = (4/5i?)f7, il being the magnitude of the angular velocity vector and R is 
the radius of the star. Transforming the electromagnetic field components in 
coordinates {x,y) and using expressions (37), (38) (39), and (40) we plot the 
trajectories of the particle for various values of the parameters M, fx, Eq, and 
Bo in figures (1) to (4); where we have taken a;(0) = 2, y{0) = 2, dx/dt \t=o= 0, 
dy/dt \t=o= and k = 1. 



7 



5 Conclusions and Summary 



In this paper we have investigated charge accretion around a slowly rotating 
compact star with a non-stationary electromagnetic field. In the region exterior 
to the star the electromagnetic field is determined by solving the generalized 
Maxwell equations (6) to (13). Then to study charge accretion in vicinity of a 
compact star, the equations of motion for a charged particle in the presence of 
the non-stationary electromagnetic field are solved numerically. For the exter- 
nal region, where the solutions are stable and convergent, our analysis shows 
that gravitational as well as electromagnetic field determine the charge particle 
trajectories. However it is the magnetic force that has the dominant effect on 
the transition from an ordered state of motion to a chaotic one. Close to the 
star particle motion is particularly sensitive to a change in the magnetic fic^ld 
of the star. Here an increase in the gravitational field strength (gravitoelec- 
tric as well as gravitomagnetic) directly causes a suppression in magnetically 
induced chaos in particle orbit (Figure 3 and 4). The process however requires 
a substantial (about 10 to 100 times) increase in the gravitoelectric field. In 
comparison a change in the gravitomagnetic parameter /i from 2 to 10 units 
introduces order in particle's orbit. On the other hand the effects of electric 
field (Figure 2) on the charged particle motion are more complicated. As the 
electric field becomes comparable to the magnetic field, the accretion trajecto- 
ries become precisely determined. Crossing this value, an increase in the electric 
field causes the charged particle to reverse its direction of motion after which 
it escapes from falling into the star. Magnetic field effects on charged particle 
trajectory displayed in figure 1 also indicate that although an increase in the 
magnetic field strength makes charge accretion chaotic, the general trend is that 
of particle in-fall due to the dominant gravitational effects. On the other hand 
in the regions close to the poles of the star, where magnetic field is particularly 
high, it is expected that chaotic trend in accretion is more dominant provided 
that the electric field effects on particle's motion are comparatively weak. 

Summarizing the above conclusions, accretion at relatively large distances 
from the surface of a compact star is more or less ordered due to the dominant 
effects of gravitational as well as electric fields. However in the star's vicinity 
accretion becomes chaotic depending on the magnetic field strength close to 
surface of the star as compared to the other forces involved, in particular the 
gravitational and the electric forces. 
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Figure Captions; 

Figure 1: Effects of the magnetic field on charge particle trajectories. The 
magnetic field strength varies from Bq = 10, 100, 1000, to 10000 units, while 
M = 1, /z = 1, and i?o = 1 in gravitational units. 

Figure 2: Effects of the electric field on charge particle trajectories. The 
electric field strength varies from Eq — 1, 10, 100, to 1000 units, while M = 1, 
/z = 1, and Bo ~ 10 in gravitational units. 

Figure 3: Effects of the gravitoelectric field on charge particle trajectories. 
The gravitoelectric field strength varies from M = 1,10,100, to 1000 units, 
while /i = 1, i?o = 1 and Bq = 10 in gravitational units. 

Figure 4: Effects of the gravitomagnetic field on charge particle trajectories. 
The gravitomagnetic field strength varies from /i = 0.05, 1, 3, to 10 units, while 
M = 1, Eo = 1 and Bq = 10 in gravitational units. 
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